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How Magnetic is the Dirac Neutrino? 
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We derive model-independent, "naturalness" upper bounds on the magnetic moments (i„ of Dirac 
neutrinos generated by physics above the scale of electroweak symmetry breaking. In the absence 
of fine-tuning of effective operator coefficients, we find that current information on neutrino mass 
implies that \\i v \ ~ 10~ 14 Bohr magnetons. This bound is several orders of magnitude stronger than 
those obtained from analyses of solar and reactor neutrino data and astrophysical observations. 
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With the current emphasis on understanding the pat- 
tern of neutrino mass and mixing and the correspond- 
ing implications for cosmology and astrophysics, it is 
also of interest to consider the electromagnetic proper- 
ties of the neutrino. The leading coupling of the neu- 
trino to the photon is the magnetic moment, ji v . The 
chiral symmetry obeyed by the massless neutrinos of 
the Standard Model (SM) requires that fj, v = 0. Now 
that we know that m„ ^= 0, however, it is interest- 
ing to ask how large one might expect the neutrino 
magnetic moment to be. In the minimally extended 
SM containing gauge-singlet right-handed neutrinos, one 
finds that \x v is non-vanishing, but unobservably small: 
/i„«3xlO- 19 /iBK/leV] Q. 

Current experimental limits are several orders of mag- 
nitude larger. Those obtained from laboratory experi- 
ments are based on analyses of the recoiling electron ki- 
netic energy T in neutrino-electron scattering. The effect 
of a non- vanishing /i„ will be recognizable only if the cor- 
responding electromagnetic cross section is comparable 
in magnitude with the well-understood weak interaction 
cross section. The magnitude of \x v which can be probed 
in this way is then given by 
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Considering realistic values of T, it would be difficult to 
reach sensitivities below ~ 10 _11 /is. The limits derived 
from studies of solar and reactor neutrinos are presently 
somewhat weaker: ^ 1.5 x 10 _10 /xb (solar) |2j and 
|/x„| < 0.9 x IO^Vb (reactor), Q. 

Limits on \x v can also be derived from bounds on un- 
observed energy loss in astrophysical objects. For suffi- 
ciently large /z„, the rate for plasmon decay into vv pairs 
would conflict with such bounds. Since plasmons can also 
decay weakly into vv pairs, the sensitivity of this probe 
is again limited by size of the weak rate, leading to Q 
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where ujp is the plasma frequency. Since (Tiuip) 2 <C m e T, 
this bound is stronger than the limit in Eq. |JT]|. Given 
the appropriate values of frujp it would be difficult to 
reach sensitivities better than lCP 12 /!^. Indeed, from 



the analysis performed in Ref. [5j , one obtains 

\fJLu\ ^ 3X 

In what follows, we show - in a general and model- 
independent way - that a magnetic moment of a Dirac 
neutrino with magnitude of the same order, or just be- 
low, current limits would be unnaturally large and would 
require the existence of fine tuning in order to prevent un- 
acccptably large contributions to m v via radiative correc- 
tions 1 . Although small Dirac neutrino masses imply very 
small Yukawa couplings, they are not inconsistent with 
observations. In order to satisfy m v < leV, we argue 
that a more natural scale for \\i v \ would be ^ 10~ 14 /is. 

Assuming that fi v is generated by some physics be- 
yond the SM at a scale A, its leading contribution to the 
neutrino mass, <5m„, scales with A as 



a A 2 \i v 

amy ~ — , 
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(3) 



where, 8m v is the contribution to a generic entry in the 
3x3 neutrino mass matrix arising from radiative cor- 
rections at one-loop order. The dependence on A 2 arises 
from the quadratic divergence appearing in the renor- 
malization of the dimension four neutrino mass oper- 
ator. Although the precise value of the coefficient on 
the right side of Eq. ® can only be obtained with the 
use of a specific model, it implies an order-of-magnitude 
bound on the size of /z„. For A ~ 1 TeV, requiring that 
8m v not be significantly larger than one eV implies that 
|m„| ^ 10~ 14 /!b. Given the quadratic dependence on A, 
this bound becomes considerably more stringent as the 
scale of new physics is increased from the scale of elec- 
troweak symmetry breaking, v ~ 250 GeV. 

The problem of reconciling a large magnetic moment 
with a small mass has been recognized in the past, and 
the quadratic dependence on A in Eq. © discussed in, 
e -g-> (a, Itl- Possible ways of overcoming this restriction 
include imposing a symmetry to enforce m„ = while 
allowing a non-zero value for [i v |8|], or employing a spin 
suppression mechanism to keep m„ small 9]. Neutrino 



1 The idea that SM-forbidden operators might contribute to m v 
through loop effects was first proposed in Ref. |(| and recently 
discussed in Ref. Q- 
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magnetic moments are reviewed in 0,^2 El; an d recent 
work can be found in (l4| . 

When A is not substantially larger than v, the con- 
tribution to 8m v from higher dimension operators can 
be important, and their renormalization due to opera- 
tors responsible for the neutrino magnetic moment can 
be computed in a model-independent way. As we discuss 
below, dimension six operators are the lowest that con- 
tribute. We shall now outline this calculation and the 
resulting constraints on Specifically, we find that 
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where 9\v is the weak mixing angle, 



/ = (l-r)-|rtan 2 w -i(l + r)tan 4 w , (5) 

and r = C_ / C+ is a ratio of effective operator coefficients 
defined at the scale A (see below) that one expects to be 
of order unity. Again taking A ~ 1 TeV, 5m v < leV, and 
setting r ~ 1, we find that \i v (for any mass eigenstate) 
should be smaller in magnitude than ~ 10 _14 //b. 

To arrive at these conclusions, we consider an effec- 
tive theory containing Dirac fermions, scalars, and gauge 
bosons that is valid below the scale A and that respects 
the SU(2) L x U(l)y symmetry of the SM. We also im- 
pose lepton number conservation. In this effective the- 
ory, the right-handed components of the neutrino have 
zero hypercharge (Y) and weak isospin. The effective 
Lagrangian involving the vr, left-handed lepton isodou- 
blet L, and Higgs doublet cj) obtained by integrating out 
physics above the scale A is given by 



E 



A n-4 



(6) 



where the n > 4 denotes the operator dimension, j runs 
over all independent operators of a given dimension, and 
fi is the renormalization scale. For simplicity, we do 
not write down the n — 4 operators appearing in the 
SM Lagrangian or the Dirac Lagrangian for the vr. At 
ri = 4, a neutrino mass would arise from the operator 
— L<f>UR, where <j> = iTii>* ■ We also omit explicit 
flavor indices on the L and vr fields. After spontaneous 
symmetry breaking (SSB) at the weak scale, 
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(7) 



so that CfO^ — > —m v VLVR with rn^ = —Cfv/\/2. Con- 
sistency with the present information on the scale of m v 
requires that \Cf \ <5x 10" 12 . 

A neutrino magnetic moment coupling would be gen- 
erated by gauge- invariant, dimension six operators that 
couple the matter fields to the SU(2)l and U(l)y gauge 



fields and By., respectively. Above the scale v, these 
operators will mix under renormalization with other n = 
6 operators that contain the L, Vr, and </> and that gen- 
erate neutrino mass terms after SSB. For this purpose, 
the basis of independent n — 6 operators that close under 
renormalization is given by 
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(fi) 



(fi) 



g-iLcfxT^VRB^ 
g 2 LT a ^ u v R W« v 
Lcj)VR 



(8) 



(4) where B^ = d„B v - d v B^ and W*, = d^W* ~ d u W£ - 
92eabcW h JV^ are the U(l)y and SU(2) L field strength 



tensors, respectively, and gi and g2 are the corresponding 
couplings. After SSB one has 



0{ 6) -> ■j*giv L o* u v R Bp, 



O 



(fi) 



g 2 v L a^v R Wl v + - 



(9) 
(10) 



Using g2 sinOw — 9i cos 9w = e, it is straightforward to 
see that the combination CfO^ + C\0^ appearing in 
£ e ff contains the magnetic moment operator 



- —vo^ vF^ 
where is the photon field strength tensor and 



-4^2 



[C!(v) + Cl(v)} 



(11) 
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Similarly, the operator O3 6 generates a contribution to 
the neutrino mass 



Sm v = -Cl(v) 



2V2A 2 



(13) 



Other n — 6 operators that one can write down are 
either related to those in Eqs. JSJ by the equations of 
motion or do not couple to F^ after SSB. It is instructive 
to consider a few illustrative examples. In particular, 
consider the following three operators: 



(fi) 



(fi) 
5 

(6) 



(14) 
(15) 
(16) 



where = + ig 2 T a W°/2 + ig^YB^/2 and where the 
sum over a = 1, 2, 3 in is implied. We may express 
Ojj in terms of O^l by first noting that 







rt^rtr 



(17) 



since Lip = = pL by the equation of motion for L. 
Then using 7^ = g^ v - ia^ we have 

Q = L^\^^R- l -L0\,L^l]<j^^ R . (18) 
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FIG. 1: Self-renormalization of O^l, denoted by the shaded 
box. Solid, dashed, and wavy lines indicate leptons, Higgs, 
and gauge bosons, respectively. 




FIG. 2: Renormalization of due to insertions of (D± 



,(6) 



FIG. 3: Self-renormalization of O 



Working out the commutator [D , D v \ in terms of B^ v 
and gives 



Y L Ol 



(6) 



o 



(6) 



(19) 



In the case of 

* 7,0* 



where Yj, = — 1 is the lepton doublet hypercharge. 

the component involving 

i?LT3. l (f> *d fi i'ii contains only the combination g 2 — 
g\B^ oc Z^ since the SM Lagrangian for the neutrino 
contains no coupling to the photon. Moreover, since O5 6 
contains a derivative acting on Vr and vr has no gauge 
interactions, it does not mix with C^ 3 under renormal- 



ization. Finally, one can show that O, 



(0) 



C3 using 



the identity t^-t^j = 2SuSjk — SijSki- Other operators that 

contain derivatives acting on the 4> may be related to O^l 
using integration by parts and the foregoing arguments. 

The one-loop renormalization of the C^ 3 is obtained 
by computing Feynman diagrams of Figs. 1-3, where 
only representative examples of the full set of graphs are 
shown. The graphs of Fig. ^ involve renormalization of 
0^2, where the shaded box indicates insertions of the 
tree-level operator. Graphs of the type shown in Fig. [21 
give renormalization of by insertions of the . At 



this order, there are no insertions of that renormalize 
the . Graphs leading to self-renormalization of O3 6 ' 
generated by gauge and X^cj)) 2 couplings are illustrated 
in Fig. [3J For the diagrams involving internal gauge bo- 
son lines, we use the background field gauge [lfjj, which 
allows us to obtain gauge-invariant results in a straight- 
forward manner. Throughout, we use dimensional regu- 
larization, working in d = 4 — 2e dimensions, and intro- 
duce the renormalization scale /j,. Due to operator mix- 
ing, the renormalized operators of! can be expressed in 



terms of the un-renormalized operators Of 1 ' via 



1/2 7 n«/2 (6) 



(20) 



1/2 1/2 

where and Z^ are wavefunction renormalization 

constants for the L and <f> respectively and where = 1 
(3) is the number of 4> fields appearing in Ofl {O^). 
In the minimal subtraction scheme that we adopt here, 
the products of renormalization constants ZJ^ Z]/ 2 Z^ 4 "' 2 
simply remove the 1 /e terms arising from the loop graphs. 

The renormalized operators are dependent on the 
scale jj. since the bare operators 



of } = zi /2 z;* /2 o 



(6) = E^i 6 i 



(21) 



must be /x-independent. The /u-dependence of the Cf(fi) 
is such that the renormalized effective Lagrangian 
does not depend on the renormalization scale. In order 
to obtain Eq. |J2J| we require the value of the Cf(fx) at 
the scale fx = v, below which the Z and M /± are inte- 
grated out of the effective theory and only the photon 
contributes to operator renormalization. Since Q v = 0, 
the latter occurs at higher order in a/47r than consid- 
ered here. The value of the Cf (v) are determined by the 
renormalization group equation (RGE) that follows from 
the requirement that C^ s be /^-independent: 



(22) 



where the anomalous dimension matrix is defined by 

d 



7fcj 



E 



(23) 



We find 



-{ax - 3a 2 ) 

8tt 2 





10 7T 



(a>i — 3a 2 ) 
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16rr 



— 6«i(ai + a 2 ) 
6a 2 («i + 3a 2 ) 
(ai + 3a 2 ) 



5 \ 



(24) 



4 



where the a, = gf/in and V(<j>) = A[(0ty) - v 2 /2} 2 . 

Using the known /3-functions that govern the [i- 
dependence of the g$ and the anomalous dimension ma- 
trix in Eq. (|2"4l we numerically solve the RGE lt2^|l for 
the Cf(n). In doing so, we find that the ^-dependence 
of the Qi has a negligible impact on the overall solu- 
tion. Neglecting the /^-dependence of the #j then al- 
lows us to obtain an analytic solution. In this ap- 
proximation we find that the combinations of constants 
C+ifi) = Cf(n)+C$(ji) andCV) = aiC?0u)-3a 2 Cf(/z) 
evolve independently. Since fi v is proportional to C+ (v) , 
the presence of a non-zero neutrino magnetic moment at 
low energy requires the physics beyond the SM to have 
generated a non- vanishing C+(A). It is then straight- 
forward to obtain Cf (//) which depends on all three of 
the Cf (A). Retaining only the leading logarithms rather 
than the full resummation provided by the RGE and 
defining C_(/x) = Cf (//) - Cf (ji) we find 



Cf(M) 



C + (A) 
C(A) 



7ln A 



(25) 



1 + 7 In 



A 



C 3 6 (A) [l- 7 33ln^ 



[C + (A) 7+ + C-(A) 7 _]ln 



A 



where 7± = (713 ±723)/2 and 7 = 3(ai + 3a2)/167r. 

Using Eqs. H12I13(I allows us to relate {i v to the cor- 
responding neutrino mass matrix element in terms of 
C±(A) and Cf (A) 



5m v — 



v 2 Cl(v) \i v 



16m e C+(v) [is 



(26) 



with C+{v) and Cf (1;) given approximately by Eqs. 
To obtain a natural upper bound on fi v , we assume first 
that Cf (A) = so that 8m v is generated entirely by 
radiative corrections involving insertions of 0^f\ . Doing 
so in Eq. (|2*rJ|) and solving for (x v /(ib leads directly to 
Eq. To arrive at a numerical estimate of this bound, 
we substitute A = 1 TeV into the logarithms appearing 
in Eq. an d obtain 



1/H < 
Mb ~ 



8 x 10 



-15 



5m u 
1 eV 



1 

171 



(27) 



It is interesting to consider the bound for the special 
case that only the magnetic moment operator is gen- 
erated at the scale A, i.e., C+(A) 7^ and C_ = 0, 
with / ~ 1. For this case, considering a nearly degener- 
ate neutrino spectrum with masses ~ 1 eV leads to the 
\Hv\ ^ 8 x lQ~ lb ^B ~ a limit that is two orders of mag- 
nitude stronger than the astrophysical bound |5j and 10 4 
stronger than those obtained from solar and reactor neu- 
trinos. For a hierarchical neutrino mass spectrum, the 
bound would be even more stringent. 



The discovery of a Dirac neutrino magnetic moment 
having a magnitude comparable to, or just below, the 
present experimental limits would imply considerable 
fine-tuning in order to maintain consistency with the 
scale of neutrino mass. Such fine-tuning could occur 
through cancellations between the C+(A), C_(A), and 
Cf (A) terms in Eq. I|25|l . While it is in principle possible 
to construct a model that displays such fine-tuning, the 
generic situation implies substantially smaller magnetic 
moments for Dirac neutrinos than are presently accessi- 
ble through observation. 



The limits one may obtain on transition magnetic mo- 
ments of Majorana neutrinos are substantially weaker 
than those for the Dirac moment. Because the transi- 
tion magnetic moment is antisymmetric in the flavor 
labels i, j, while the mass matrix is symmetric, 5m v must 
be higher order in or involve insertions of the Yukawa 
couplings. 
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